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Extracting the intrinsic switching field distribution in perpendicular media: a 

comparative analysis 

Michael WinklhofeiE and Gergely T. Zimanyi 
Department of Physics and Astronomy, UC Davis, One Shields Avenue, Davis, CA 95616, USA 

We introduce a new method based on the first-order-reversal-curve (FORC) diagram to extract the 
intrinsic (microscopic) switching-field distribution (SFD) of perpendicular recording media (PRM). 
To demonstrate the viability of the method, we micromagnetically simulated FORCs for PRM with 
known SFD and compare the extracted SFD with the SFD obtained by means of two different 
methods that are based on recoil loops, too, which however rely on mean-field approximations 
and assumptions on the shape of the SFD. The FORC method turns out to be the most accurate 
algorithm over a broad range of dipolar interaction strengths, where the other methods overestimate 
the width of the SFD. 

PACS numbers: 85.70.Ay,75.50.Ss,75.50.Tt,75.60.-d 



I. INTRODUCTION 

The quality of recording media depends crucially on 
the intrinsic (microscopic) switching-field distribution 
(SFD) of the media particles, which determines both 
magnetic stability and attainable recording density. It is 
straightforward to obtain the SFD of a diluted magnetic 
system by taking the derivative of the DC demagnetiza- 
tion (DCD) curve. In the case of high-density magnetic 
recording media, where magnetic interactions between 
the media particles are not negligible, the problem of ex- 
tracting the SFD from bulk magnetization curves cannot 
be solved rigorously any more and the shape (and to a 
lesser extent the location) of the extracted SFD will de- 
pend on certain model assumptions. 

Two conceptually different algorithms 1 ^ have been 
suggested recently to extract the intrinsic D(H S ) of per- 
pendicular recording media (PRM) from macroscopic 
magnetization curves using a set of recoil loops. The 
analysis technique by van de Veerdonk et ali assumes 
a constant effective demagnetization factor N e g for des- 
hearing recoil loops, from which the DCD curve is ex- 
tracted. Extraction of the DCD curve and deshearing 
the recoil loops are performed simultaneously in order to 
arrive at a self-consistent solution!. Although the algo- 
rithm converges after a couple of iterations, it is not clear 
how reliable the self-consistent solution may be under the 
assumption of an effective demagnetization factor inde- 
pendent of the magnetization. Another drawback of the 
method is that it requires assumptions about the shape 
of the SFD, which of course is not known. 

The second algorithm, referred to as AH(M, AM)- 
method 2 , is a generalization of the AH method orig- 
inally proposed by Tagawa and NakamuraA The 
AH(M, AAf)-method overcomes the restriction of a con- 
stant value of N e g, assuming an effective field of the form 
H eS (H a ,M) = H a + Hi(M), where H a and Hi denote 
the applied and internal field, respectively, and M is the 
magnetization in the direction of H a . The AH(M, AM)- 
method approximates interactions on the mean-field level 



in a sense that all microscopic magnetization configu- 
rations M(r) representing the same macroscopic mag- 
netization value M produce the same volume-averaged 
internal field H{(M). An implicit assumption underly- 
ing the AH(M, AM)-method is that each particle acts 
a square hysteron. Then, the lower branch of the major 
loop M~(H a ) can be represented as: 



M-(H a ) = -M sat + 2 



H etf (Ka,M) 



D{H s )dH s , (1) 



where D(H B ) is the SFD. Any recoil loop M J (H r , H a > 
H r ) originating from the upper branch M + at H r can be 
written as 



M J (H ri H a ) = M + (H r ) + 2 



H cll (H a ,M) 



D(H s )dH s . (2) 



In order for Eq. [5] to hold, two conditions must be 
met: firstly, M J (H r , H a > H r ) has to saturate at fields 
I -Ha | < \H r \ and secondly, recoil loops must not cross 
each other. Then the coercivity distribution can be as- 
sumed to consist of disjunct segments D(H s )dH s (non- 
interacting hysterons). The inverse i -1 of the cumulative 
distribution 



1(H) 



H 



D(H s )dH s 



(3) 



for a given magnetization value M can then be obtained 
by taking the difference in field position, AH(M), be- 
tween M~(H a ) and recoil loop M J (H r ,H a > H r ), 



AH(M, AM) = H a (M~ 



= r 



M 



- I 



H a {M J ) = 
M sat 



(4) 



(M - AM) 



where AM = M + (H I ) - (-M sat ). By fitting the 
AH(M, AM) curves against the inverse of a certain pa- 
rameterized distribution function, the key features of the 



2 



SFD can be extracted^. In case I is a normal distribution 
centered at H c , Eq. 0] writes to 



AH(M, AM) 



V2a [erf" 1 (to) - erf _1 (m - Am)] , 



(5) 

where lower case to's are the magnetization values rela- 
tive to the saturation magnetization M saA - 

The third method in our comparative analysis is based 
on the first-order-reversal curve (FORC) distribution 4 ^, 
defined as 



ld 2 M(H T ,H a 
2 dH,dH a 



(6) 



Although there is no principal difference between a 
FORC and a recoil loop, the numerical determination 4 
of p from a set of FORCs requires recoil loops measured 
on a grid equidistant in H r and H a . 

For an assemblage of square hysterons, the FORC dis- 
tribution is identical to their Preisach distribution and 
we can obtain their SFD as 



/oo 
p(H c ,H h )dH h 
-OG 



(7) 



where H c = (H a - H T )/2 and H h = (H a + H t )/2 are the 
coordinates of the Preisach plane, defined by the micro- 
scopic switching field H c and the bias field H\>. As with 
the AH(M, AM) method, the FORC method is reliable 
as long as the media particles can be reasonably well de- 
scribed by square hysterons. As opposed to the two algo- 
rithms above, however, the method based on Eq.0 relies 
neither on a mean-field approximation nor on a certain 
model form for D(H S ) and therefore imposes the fewest 
constraints on the data. 



II. COMPUTATIONAL DETAILS 

The three methods presented above are best tested on 
a system with known SFD. For this purpose, we mi- 
cromagnetically computed sets of FORCs for a PRM, 



using the OOMMF code 6 



l.lb2). The simulated 



medium typically consisted of N ~ 2 • 10 3 particles, 
arranged on a regular grid (mesh size s = 5 nm), in- 
scribed in a circle of 250 nm diameter. We chose a cir- 
cular boundary to minimize the effects of corners. The 
easy axes of the particles are all oriented the same way, 
roughly perpendicular to the surface (9k = 89.42 deg) 
so as to avoid numerical problems that may arise when 
the applied magnetic field H (Oh — 90 deg) is exactly 
collinear with the easy direction. We modified the mod- 
ule maganis . cc to produce normal distributed values of 
the uniaxial magnetocrystalline anisotropy constant, K u 
such that (t Hk = 0.1 (H K ), with H K = 2 K u /p Q M sat . ft 
has been showni that a small amount of intergranular ex- 
change helps to reduce the increase in the recording tran- 
sition parameter due to a small distribution in Hk- We 
therefore set the nondimensional intergranular exchange 
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FIG. 1: Comparison of the three methods in terms of the ex- 
tracted a* of the SFD as a function of the dipolar interaction 
strength (f>q. The relative strength of the exchange coupling 
decreases from left to right (/i cx = 0.046 for all values of 0q). 
The intrinsic SFD is a Gaussian distribution with a = 0.1 
(thin horizontal line). The coercive force is represented by 
the grey line without plot symbols (_ff c scale on the right). 



coupling constant /i C x = A/(K u s 2 ) to /i cx = 0.046 in 
all our simulations. In order to systematically explore 
the effects of dipolar interactions, we vary the dipolar 
interaction strength by using a scaling factor <f> < 1 for 
-Msat, where <p represents the volume fraction of mag- 
netic material in each cell. The uniaxial magnetocrys- 
talline anisotropy constant, K u , is scaled simultaneously 
to keep the value of the microscopic coercivity Hk con- 
stant as 4> varies. We normalize all magnetic fields by Hk 
and use a mean-field type equation to define the material- 
independent dipolar interaction strength as 4>q = (ft/Q, 
i.e., 



M 



a- 



= h a + at 



(8) 



QM sat 

where a is the mean-field parameter and Q denotes the 
quality factor Q = 2K u /(^i M 2 at ). For hep Cobalt at 
room temperature, Q = 0.422. 



III. RESULTS AND DISCUSSION 

Figure ^ shows the estimated width a* of the SFD as 
extracted with the three algorithms tested here. In the 
hypothetical limit case of a PRM controlled solely by 
exchange-coupling (4>q < 0), the AH(M, AM) method 
is able to capture the intrinsic width of the SFD with 
great accuracy. For PRM largely dominated by dipo- 
lar interactions, however, the AH(M, AM) method is 
consistently less accurate than the FORC method, which 
renders the least deviation of a* from a for </>q > 0.15. In 
contrast, the a* extracted with the N e g method always 
show the largest deviations from a. Importantly, the N e g 
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FIG. 2: (Color online). FORC diagram p(H Cl H b ) (above) 
obtained for a PRM with <j>Q = 0.35, h eyi = 0.046 and the 
corresponding SFD extracted on the basis of Eq. (below). 
The width a* = 0.107 is determined from fitting the ex- 
tracted SFD (solid line) with a normal distribution centered 
at H c — 0.76 Hk (dashed). The derivative of the DCD curve 
(dash-dotted) has a distinct asymmetry. The intrinsic SFD 
(dotted) with a = 0.1 is shown for comparison. All curves are 
normalized to unit area. With the AH{M, AM) method, a* 
was obtained as 0.133 (original recoil loops and AH(M, AM) 
curves available online). 



method does not yield self-consistent solutions any more 
for 4>q > 0.25, where the asymmetry of the DCD curve 
is too pronounced to be properly described by a normal 
distribution (Fig. More consistent solutions in this 
regime can be found with a log-normal SFD cut-off at 
high fields. Caution should therefore be taken to not in- 
terpret asymmetric DCD curves prematurely as evidence 
of an asymmetric intrinsic SFD. The FORC method on 



the other hand captures the intrinsic shape of the SFD to 
a very good degree despite the shift in location (Fig. |5J . 

All the three methods are primarily concerned with 
finding the right scale parameter (i.e., a* or FWHM) 
of the SFD distribution, while the location parameter is 
considered to be invariant. As can be seen in Figure ^ 
the observed values of H c are progressively shifted to 
lower values with increasing <pQ, which is due to dipolar 
interactions between the media particles, deflecting the 
effective field away from the easy axes. According to the 
Stoner-Wohlfarth relationship, 

H B (il>) = H K ((cos^) 2 / 3 + (sin^) 2 / 3 )~ 3/2 (9) 

observed coercive forces of ~ 0.8Hk suggest that a large 
fraction of grains experience a local effective field that 
deviates by some 5% from the applied field direction, 
giving rise to curvilinear hysterons, in other words, re- 
versible magnetization processes. In the FORC dia- 
gram (Fig. [2J , curvilinear hysterons manifest themselves 
in the form of the small negative region centered at 
(H c ,Hb) — (1, — 0.05)i?K- This way, one can determine 
directly from the FORC diagram if and to what degree 
the model assumptions are not strictly met. 

Since all methods tested here rely on square hysterons, 
their underlying assumption starts breaking down with 
increasing dipolar coupling <f>Q, albeit at different values 
of <j>Q. More importantly, the fact that the methods do 
not all start to fail at the same point also shows that the 
presence of reversible magnetization processes is not the 
most crucial limitation to a method's applicability, which 
rather is restricted by imposing constraints on the inter- 
nal field distribution. Compared to the N e g method, the 
AH(M, AM) method puts a less tight constraint on the 
expected magnetization curves and so has a larger range 
of applicability. The very absence of any such approxima- 
tions makes the method based on the FORC distribution 
the most robust algorithm and therefore most suitable 
for characterizing ultra- high-density PRM. 
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Micromagnetically computed FORCs for <j)Q = 0.35. 
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Average deviation of expected AH(M, AM) curves (ac- 
cording to Eq. |3J from AH{M, AM) dataset as a 
function of a. The best fit is obtained for ctfjt = 0.133, 
which defines a* . 
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Comparison between expected AH(M, AM) curves 
(solid lines) for a* = 0.133 and AH(M, AM) dataset 
(dots). Not all AH(M, AM) curves are shown to keep 
the picture tidy. 



